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Abstract 

By considering analytical expressions for the self-energies of intervalley and intravalley phonons 
in graphene, we describe the behavior of D, 2D, and D' Raman bands with changes in doping (fi) 
and light excitation energy (El). Comparing the self-energy with the observed [i dependence of the 
2D bandwidth, we estimate the wavevector q of the constituent intervalley phonon at fivq ~ El/1-6 
(v is electron's Fermi velocity) and conclude that the self-energy makes a major contribution (60%) 
to the dispersive behavior of the D and 2D bands. The estimation of q is based on a concept of 
shifted Dirac cones in which the resonance decay of a phonon satisfying q > oj/v (u is the phonon 
frequency) into an electron-hole pair is suppressed when [i < (hvq — hu)/2. We highlight the 
fact that the decay of an intervalley (and intravalley longitudinal optical) phonon with q = lo/v is 
strongly suppressed by electron-phonon coupling at an arbitrary [i. This feature is in contrast to 
the divergent behavior of an intravalley transverse optical phonon, which bears a close similarity 
to the polarization function relevant to plasmons. 
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The Raman spectrum of graphene has two prominent peaks called the G and 2D (or 
G') bands that are very informative characterization tools. The 2D band at ~ 2600 cm -1 
has been used to distinguish a single layer from graphene layers.- The G band at ~ 1580 
cm -1 can be used to determine whether or not the position of the Fermi energy ft is close 
to the Dirac point, since the width broadens when fx ~ 0.-"- By contrast, the 2D bandwidth 
sharpens when fx ~ 0.— ^ What is the origin of the difference between the fx dependencies of 
the G and 2D bands? 

As illustrated in Fig. [U^a), the presence (absence) of a resonant process by which the 
phonon decays into a real electron-hole pair, enhances (suppresses) the spectral broadening. 
Because the G band consists of T point phonons, a direct transition is a unique decay 
channel that conserves momentum. Thus, the fx dependence of the G bandwidth is readily 
understood in terms of the Pauli exclusion principle.- - - Meanwhile, the 2D band involves two 
near K point (intervalley) phonons, 8,9 and the spectral broadening is induced by an indirect 
transition that crosses two valleys, as shown in Fig. QJb). The presence or absence of a 
resonance decay channel for a phonon with a nonzero wavevector is the key to answering 
the question posed above. In this paper, we provide the answer in a unified manner by 
translating the Dirac cone. 

Figure [T](b) shows that an intervalley phonon (zigzag line) can change into an electron-hole 
pair (loop) as a result of an electron-phonon interaction. The wavevector of an intervalley 
phonon is written as kp + q, where kp is a wavevector pointing from the K point to the K' 
point and q (= |q|) is much smaller than |kp|. Suppose that a hole is located at k measured 
from the K point, then the wavevector of the electron is given by k + (kp + q) because of 
momentum conservation. As a result, the wavevector of the electron measured from the K 7 
point is k + q, and the energies of the hole and electron are given by —hvk (= — 7u>|k|) and 
fa;|k + q|, respectively, where v (~ 10 6 m/s) denotes the Fermi velocity. Hereafter, we use 
units in which Ti — 1. 

An indirect transition between two valleys can be regarded as a "direct" transition by 
translating the Dirac cone at the K' point to — (kp + q) as shown in Fig. Hfb). With 
the shifted Dirac cones, it is easy to capture the essential feature of the broadening of a 
phonon. When fx = 0, we see in Fig. [T(c) that there is an energy gap, vq, between the 
conduction and valence bands. This energy gap precludes a phonon with frequency to (< vq) 
from decaying into a real electron-hole pair. On the other hand, when sufficient doping is 
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FIG. 1: (a) The /i dependence of G band broadening (frequency u)q). (b) An electron-hole pair 
between two valleys. The Dirac cone at the K (K 7 ) point is indicated in black (red). With the 
migration of the Dirac cone, the electron-hole pair creation process is viewed as a direct transition, 
(c) The /i dependence of the broadening of an intervalley phonon is different from that of the G 
band. Note that spectral broadening of an intravalley phonon can be discussed with the replacement 
K' -> K (or K -> K'). 

achieved as shown in Fig. Q^c), the phonon can decay into an intraband electron- hole pair. 
This intraband decay channel results in spectral broadening. When q = 0, the two Dirac 
cones are merged into one [inset in Fig. QJb)] and the energy gap vanishes. Then, it is clear 
that the broadening of the q = phonon bears similarities to that of the G band." The 
/i dependence of the broadening of an intervalley phonon with vq > u differs greatly from 
that of the G band, and the concept of the shifted Dirac cones is useful for understanding 
the fi and q dependencies of the broadening in a unified manner. 

More detailed information about the broadening can be obtained by calculating the self- 
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energy. The self-energy of the intervalley phonon with q and u (> 0) is defined by 

n,(g, = 9cp — £ L a; + W fc-Hk + q|+»6 I 1 ~ SS "TkTqF ) 

s,s' k 11/ 

In Eq. (HJ, s (= ±1) and s' (= ±1) are band indices, = lim^^l + e ^(«W-M))-i i s 
the Fermi distribution function defined at zero temperature and with a finite doping /i, and 
e is a positive infinitesimal. We can assume fi > without losing generality because of 
particle-hole symmetry. The factor g ep denotes the electron-phonon coupling strength, (p 
denotes the polar angle between k and q, and the term, g cp x (1 — ss ^l^"^ ), is the square 
of the electron-phonon matrix element for the intervalley phonon, which will be discussed 
later. The broadening and modified frequency are given by — ImILj(g, u) and cu + ReU^(q, u>), 
respectively. 

In the continuum limit of k, the broadening normalized by g ep leads t^~ 
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where 6 X denotes the step function satisfying 9 X > = 1 and 9 X<0 = 0, and g(x) = log(x + 



\Jx 2 — 1). Figure |5fa) shows a 3-dimensional (3d) plot of — ImILj(g,u;) as a function of vq 
and /i when u = 0.2 (eV), which corresponds to the Debye frequency of carbon (wo). As 
indicated in Fig. E|a), a line node appears for vq = ujd- In Eq. (J2J), this line node is a 
critical line separating the two terms, which are proportional to 9 W - Vq and 9 vq ^, and it 
can be shown that the first (second) term originates from the contributions of interband 
(intraband) electron-hole pair creation. For example, the q = phonon satisfies 9 vg - u = 0, 
and only the interband electron-hole pairs cause spectral broadening. The first term in 
Eq. ([2j) leads to — ImII jtt<u; /2(0, u) = ir 2 uj and ImILj >W/ / 2 (0, oS) = 0, which are consistent with 
the behavior of the G band.-~- Contrastingly, for the phonon satisfying vq > u, we can 
confirm that from the second term of Eq. (J2J) spectral broadening is possible only when 
there is sufficient doping, namely when \i > ^2^L j s satisfied. A sharp step appears at 
v q — 2fM + uj d , as indicated in Fig. Eta). In Fig. [2(b), we plot — ImlL^g, ud) as a function 
of \i to show more clearly the q dependence of the broadening. It is seen that for vq > 0.2 
[red curves], — ImII M (g, cod) is suppressed when the Fermi energy is close to the Dirac point 
{jj, < vq Z " 2 ), and broadening appears when ji > vq ~ ' 2 . 



Interband 




FIG. 2: (a) 3d plot of — ImU^q, cjd). The variables vq and \x are given in units of eV. The cross 
section of (a) for different vq values (b), and for different \i values (c). In (c), vq is proportional to 
the light excitation energy El. 

Because the Raman 2D band consists of two intervalley phonons satisfying vq > u, 9 the 
suppressed broadening when /x ~ also holds for the 2D band. Das et al. A have shown 
that the 2D bandwidth sharpens when \i ~ 0. A suppressed broadening of the 2D band 
(2u q ~ 0.32 eV) has also been observed when /i < 0.4 eV in a recent experiment reported 
by Chen et a/.,- from which we estimate the q value to be vq ~ 0.96 eV using 0.4 ~ 16 . 
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The validity of this estimation (vq ~ 0.96 eV) can be further investigated by changing q. 
In Fig. Efc), we show the plot of — Imil M (g, ud) as a function of q for different \x values. For 
vq ~ 0.96, increasing q would cause the broadening to decrease (increase) when \i = 0.4 (0.8) 
eV. Because vq is related to light excitation energy El through momentum conservation,- 
the broadening can also depend on E L - If we assume E L = avq, a ~ 1.6 is obtained as a 
fitting parameter because El = 1.58 eV is used in the experiment.- A similar a parameter 
value (a ~ 1.3) can be obtained by calculation.— 

A 3d plot of the real part of the self-energy, ReILj(g, cod), is shown in Fig. [3]^a). The plot 
is based on the analytical expression of Kell^q, u) given by^ 

ReILj(g, oj) = Anji 
+ n^/cu 2 — v 2 q 2 9 l 

+ -k\/ v 2 q 2 — uj 2 9 vq - u ^9v q -u> 

(3) 

For the q = phonon, Fig. EJ^a) shows that the softening is maximum at fi = 0.1 (eV). 
Equation [3] is simplified in the limit of q — > 0, as 

ReILj(0, u) ~ Ana + iru log 

and the large softening is caused by the logarithmic singularity at fi = ud/2. This feature 
is exactly the same as the Kohn anomaly 10 of the G band.-&^ Figure |3]^a) shows that the 
logarithmic singularity is removed gradually as we increase q from zero. [The logarithmic 
singularity is obscured by charge inhomogeneity^] It also shows that when fi is sufficiently 
large the real part increases linearly with [i as RelL^g, ojd) ~ Anfi for an arbitrary q value. 
Interestingly, ReILj(g, w^) increases as we increase q (or El), even for a fixed /i value. This 
feature is more clearly seen in Fig. M^c), and suggests that the self-energy contributes to 
the dispersive behavior of the 2D (or D) band:— ~— the 2D band frequency increases linearly 
with El {doj2T>/ 9El — 100 cm _1 /eV).— ^ If we use g ep = 5 cm -1 , which is obtained from 
the broadening data published by Chen et al.,— the self-energy can account for ~60% of the 
dispersion because 2 x ReU^ (q,u q ) ~ 2g cp Tr 2 vq = 2g ep n 2 E L / a and 2g ep n 2 /a = 61.6 cm -1 . 

In Fig. |3|a), for a fixed vq that is larger than lod, ReII At (g, ujd) undergoes two disconti- 
nuities at [i = vq ~^ D and vq+ ^ D . A modest softening appears as /i approaches vq ~^ D from 
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FIG. 3: (a) A 3d plot of ReILj(q, ujd)- The variables vq and /j, are given in eV. (b) The cross 
section of (a) for different vq values, and (c) for different \i values. Note that Rell^g, uj) does not 
include the q dependence of the bare frequency. 

zero. This is consistent with the observations by Das et a/.—, Chen et al.—, and Araujo et 
alrl- showing that the 2D band frequency remains almost constant (disregarding a small 
modulation of about 8 cm -1 ) when the Fermi energy is near the Dirac point. On the other 
hand, the 2D band frequency exhibits a slight hardening of ~ 2 cm -1 in the observation 
reported by Yan et al- We consider that the data actually show that the 2D band frequency 
does not depend on doping because the observed small amount of hardening is within the 
spectral resolution (2 cm -1 ). 2 

As we increase \i further, ReIT At (g, ujb) undergoes slight hardening and subsequent soft- 
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ening until vq+ ^ D , These features can also be seen in Fig. EJ^b). The discontinuities of 
KeHn(q, ui) can be explained by the perturbation theory: the energy correction by a virtual 
state with energy e is proportional to 

1 



co — e 

-1 ir 



(5) 



Because the sign of (u — e) is positive (negative) when e < uo [e > u), the lower (higher) 
energy electron- hole pair makes a positive (negative) contribution to ReIl A1 (g, u).— Therefore, 
when vq > u, softening is induced by the doped carriers since the energy of a virtual state 
is approximately given by vq, which is larger than u, and thus l/(u> — vq) < is satisfied 
[see Fig.|H(a)]. In fact, the energy e corresponds to v\k + q| — vk in Eq. ([I]) and e ~ v q when 
k ~ 0. The softening magnitude is tiny as shown in Fig. [3]Ja) and (b) because the electron 
density vanishes at the Dirac point. When = vq ~ UJ , an intraband electron-hole pair with 
e < oj can start to be excited [see Fig. H](b)], and this electron-hole pair causes hardening. 
Note that some of the doped carriers satisfying < ^ < uq±u_ con t r ibute to the softening, 
and the hardening is partly cancelled by the softening. The details of the cancellation are 
determined by the ip dependence of the electron- phonon coupling term, 1 — ss ' ^^^^ > m 
Eq. (CQ). Because the intraband transition satisfies ss' = 1, the matrix element vanishes when 
tp = and thus hardening dominates softening (unless u is negligible compared with vq). 
When /i approaches vq ^ u , the Pauli exclusion principle forbids the occurrence of some of the 
intraband transitions that contribute to the hardening, which accounts for the appearance 
of the softening. For /i > ^j^, the frequency exhibits hardening due to the suppression of 
the softening induced by interband (ss' = — 1) virtual electron- hole pairs. 

As we have seen, a phonon's self-energy can be very sensitive to the electron-phonon 
matrix element. In particular, the sign of the coefficient of ss' in 

, k + qcostp 

is critical in determining the behavior of the self-energy. In fact, if the minus sign is replaced 
with a plus sign as follows 

,A; + gcos^ 

the corresponding self-energy exhibits a singularity at vq = u>. In addition, the real part 
exhibits softening as ~ —Anfi when vq > u because the matrix element is maximum when 
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FIG. 4: The projection of the Dirac cones, (a) Slight doping (red) causes a softening, while (b) a 
heavy doping (blue) causes hardening as well as softening. 

if = and thus softening dominates hardening. These features are in marked contrast to the 
fact that the imaginary part of the self-energy given by Eq. fl6]) exhibits a nodal structure 
at vq — oj and that the real part exhibits hardening as ~ 47r/i. The self-energy with Eq. ([7]) 
corresponds to the self-energy of the Coulomb potential known as the polarization function, 
and the singularity at vq = u is important for plasmons in graphene.— ^ 

Intravalley longitudinal optical (LO) and transverse optical (TO) phonons are related to 
the intervalley phonon and plasmon. The corresponding elements of the electron-phonon 
interactions are given by2£ 

,k + qcostp .2/csm 2 v? 
~ ~\k + qT + TkTqf' <L0) ' <8) 
lW * + ,cos<f _ 33 , 2fcsinV 
|k + q| |k+q| 

The first two terms for the LO [TO] phonon are the same as Eq. (jSJ) [Eq. ([7j)]. By constructing 
an analytical expression for the self-energy of the LO and TO phonons, 27 we visualize the 
self-energies in Fig. [5j In Fig. [5(a), we show the imaginary part of the LO phonon. A 
notable feature of Fig. [3(a) is that for vq > u>d, the broadening increases as /x is increased. 
This is a sharp contrast to the broadening of the intervalley phonon, which is suppressed for 
heavy doping (see Fig. 12(a))- For the real part shown in Fig. [5(b), a discontinuous feature 
caused by the last term in Eq. (jSJ) is clearly seen at vq = cup. Interestingly, the TO phonon 
has some similarities to the polarization function: the existence of a singularity and the 
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frequency softening for vq > ojd, as shown in Fig. [5]^c) and (d). It is instructive to compare 
the self-energy of an intravalley LO phonon with that of a TO phonon. The self-energy 
of the LO phonon is the same as that of the TO phonon for the T point q = 0, and their 
difference is highlighted for nonzero q values (q > u/v). The difference between the LO 
and TO phonons will be useful in allowing us to determine the optical phonon (LO or TO) 
composing the D' band. 




FIG. 5: 3d plot of U^(q,C0D) for an intravalley LO phonon (a,b) and TO phonon (c,d). (a,c) the 
imaginary part, and (b,d) the real part. The variables vq and \i are given in eV. 

Coulomb interactions among electrons which we did not consider in this paper might make 
an effect on the self-energies of phonons.— Attaccalite et al. point out the importance of 
vertex-corrections to electron-phonon coupling by electron-electron interaction.— However, 
their results concern the phonon at the exact K point, which corresponds to the q = 
phonon, and the q = phonon has nothing to do with the experimentally observed 2D 
band phonon. Our main results remain unchanged even if we include such a correction, 
because our approach - shifting the Dirac cones - is based only on momentum conservation 
and does not depend on any dynamical details. Moreover, since we have determined the 
electron-phonon coupling using experimental data, such effects, if any, are all included. 
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In conclusion, employing a concept of shifted Dirac cones, we clarified that a phonon 
satisfying vq > u does not decay into an electron-hole pair when fi < (vq — u)/2. This is 
a general consequence that is independent of the details of electron-phonon coupling and 
that can be applied to both inter and intravalley phonons. Based on the self-energy, which 
includes the effect of electron-phonon coupling, we estimated the q value of the 2D band at 
vq ~ 1 (eV) by referring to recent experimental data on the \x dependence of broadening. 
This value vq ~ 1 (eV) also suggests that about 60% of the dispersive behavior can be 
attributed to the self-energy. Since vq is proportional to El, the q dependence of the self- 
energy may be explored by using a tunable laser, without changing \i by controlling the gate 
voltage. For example, the Fermi energy position of graphene can be determined from the 
El dependence of the broadening. Several anomalous features have been pointed out in the 
self-energies for intravalley LO and TO phonons. The differences between the LO and TO 
phonons will be useful for specifying the mode and q value of the D' band. 
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